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ABSTRACT 

We present the results of a parsec-scale polarization study of three FRI radio 
galaxies - 3C66B, 3C78 and 3C264 - obtained with the Very Long Baseline 
Array at 5, 8 and 15 GHz. Parsec-scale polarization has been detected in a large 
number of beamed radio-loud active galactic nuclei, but in only a handful of the 
relatively unbeamed radio galaxies. We report here the detection of parsec-scale 
polarization at one or more frequencies in all three FRI galaxies studied. We 
detect Faraday rotation measures of the order of a few hundred rad m~^ in the 
nuclear jet regions of 3C78 and 3C264. In 3C66B polarization was detected at 
8 GHz only. A transverse rotation measure gradient is observed across the jet 
of 3C78. The inner-jet magnetic field, corrected for Faraday rotation, is found 
to be aligned along the jet in both 3C78 and 3C264, although the field becomes 
orthogonal further from the core in 3C78. The RM values in 3C78 and 3C264 
are similar to those previously observed in nearby radio galaxies. The transverse 
i?M gradient in 3C78, the increase in the degree of polarization at the jet edge. 
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the large rotation in the polarization angles due to Faraday rotation and the 
low depolarization between frequencies, suggests that a layer surrounding the 
jet with a sufficient number of thermal electrons and threaded by a toroidal or 
helical magnetic field is a good candidate for the Faraday rotating medium. This 
suggestion is tentatively supported by Hubble Space Telescope optical polarimetry 
but needs to be examined in a greater number of sources. 

Subject headings: galaxies: individual (3C66B, 3C78, 3C264) — polarization 



Introduction 



Knowledge of the ambient medium near the jet launching site in the central parsecs of an 
active galactic nucleus (AGN) can shed light on why these AGNs produce j ets with different 
morphologies: low-luminosity Fanar off- Riley type I radio galaxies (FRI; iFanaroff fc Riley 
19741 ) have jets that flare into diffuse radio lobes on scales of a few hundred parsec, while 
the narrow, coUimated jets in the high-luminosity FR type II radio galaxies terminate in 
bright hot spots and create radio lobes via plasma backflow. Although the origin of the FR 
dichotomy is an open question, and a number of differences in the two classes have been cited 
as potential causes, differences in the nuclear environments of FRIs and FRIIs, which in turn 
influence the evolu tion of the jet due to interaction and entrainment , are widely believed to 
be im portant (e.g., 
1996h . 



Begelmanlll982l : iBaum et al.lll995l : lBicknelllll995l : lOe Young|ll996l : iLaing 



Very Long Baseline Interferometry (VLBI) is currently the only technique by which 
the parsec-scale environments of radio galaxies can be probed. Multi-frequency polarization 
sensitive VLBI experiments can determine the Faraday rotation measure {RM) arising in 
the vicinity of the parsec-scale radio jets. Radiation propagating through a magneto-ionic 
medium has its plane of polarization rotated from its intrinsic value xo to x following the 
relation 



where 



X = Xo + RMX' 



RM = 812 / n,B»dl rad m^ 



(2) 



rig being the electron density in cm~^, the net line of sight magnetic field in milli- 
Gauss, and dl an element of path length through the plasma in parsecs (see review by 
Gardner &: WhiteoakI 119661 ). The net parsec-scale i?M is dominated by contributions from 
the nuclear environment of the AGN in its host galaxy and our own Galaxy; the Galactic 
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contribution is usually estim ated using arcsecond-res olution observations at multiple wave- 



lengths near 18-20 cm (e.g., iRudnick fc Joned Il983l ). VLBI rotation measures have been 



observed and studied in a number of relativistically beamed AGNs, i.e., radio-loud quasars 
and BL Lac objects. Such studies have in dicated that contributions to the RM from regions 
close to the radio jets are dominant (e.g., IZavala fc Tayloiil2003l . |2004| ). 



Radio galaxies are beheved to be the plane-of-sky counterparts of quasars and BL Lacs 
(lUrry &: Padovanilll995[ ). Since their cores and jets are not as strongly Doppler boosted and 
are therefore radio faint, there have been relatively few VLBI polarization stu dies done on 
thera, and polarization has been detected on parsec-scales in onlv a handful (jTaylor et al. 
2OO1I : iMiddelberejbooi IZavala fc Tavlojbooi iKharb et allboosl : IXavlor et allbooeh . Parsec- 
scale rotation measures have only been estimated for the FRII radio galaxies 3C111, 3C120 
and 3C166, and the FRI radio galaxies M87 and 3C84. Here we augment these data with 
the parsec-scale polarization images of three FRI radio galaxies — 3C66B, 3C78 and 3C264, 
along with the RM images of the latter two. 

We have adopted a cosmology in which Hq=71 km s~^ Mpc~\ ^^=0.27 and Q\=0.73. 
The spectral indices a are defined such that the flux density S^, at frequency ly is oc z/~". 



2. The FRI galaxies under study 



We observed four FRI radio galaxies from a larger sample of radio-loud sourc es (IKharb fc Shastri 



2004 ) that possessed optical nuclei in th eir Hubble Space Telescope (HST) images (IChiaberge et al. 



19991 . 12OOOI : IVerdoes Kleiin et al.ll2002l ). and had a correlated VLBI flux density >100 mJy, 



with global V LBI at 8.4 GHz in March 2002. Parsec-scale polarized emission was detected 
in all of them (IKharb et al.ll2005l ). New observations of three of the four FRIs that showed 
polarization in their parsec-scale jets, viz., 3C66B, 3C78 and 3C264, were obtained at 5.0, 
8.4 and 15.3 GHz in September 2005 with a global VLBI array including the Very Long 
Baseline Array (VLBA) and the Effelsberg antenna, and are the focus of the present work. 

3C66B resides in an elliptical galaxy at a distance of 91 Mpc that forms a part of a 
dumbbell pair close to the Abell cluster 3 47. It is one of the f ew FRI ra dio galaxies in whic h 
the dominant radio jet ha s an infrared (ITanslev et al.ll2000l ). optical (jFraix-Burnetl 119971 ) . 
and an X-ray counterpart (IHardcastle et al.ll2001l). Kiloparsec-scale o bservations of 3C66B 



suggest the jet inclination to be less than ~ 53° (IHardcastle et a. 
for the VLBI jet of ~ 45° has been inferred by iGiovannini et al. 



I99J, while an inclination 



(I2OOII ). X-ray observations 



have revealed the presence of copious amounts of hot gas, with which the kiloparsec radio 
jet seems to be interacting (ICroston et al.ll2003l ). 
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3C78 resides in the isolated, small E/SO galaxy NGC1218 flMatthews et al.lll964l ) at a 
distance of 124 Mpc. HST o bservations have re v ealed that the opti cal galaxy has a face-on 
dust disk and an optical jet (ISparks et al.lll995l ). iBaum et al.l (Il988l ) have reported that the 
innermost emission-line gas is elongated roughly in the direction of t he radio jet, following 
the "alignment effect" fe.g.. iMcCarthv et al.l 1 19871 : |Pri von et aliboosh . 



The radio morphology of 3C264 is of the narrow- angle-tail type ( lO'Dea fc Owerull985l ). 
It is hosted by an SO galaxy 93 Mpc away that lies in a dense part of the Abell cluster 
1367 ( SchombertI 19861). Like 3C66B and 3C78, the radio jet in 3C264 also has an optical 
counterpart ( Cranel 1993). iJ6T i mages show that the optical jet is projected onto a face- 
on dust dis k (Ide Koff et al.ll2000l). A jet inclination angle of ~ 50° h as been inferred on 



kiloparsec- (IBaum et al 



19971 . Il998l ) and parsec-scales (iGiovannini et al.ll200ll ). 



3. Observations and Data Reduction 

Observations were carried out with the ten elements of the VLBA and the 100-m Ef- 
felsberg telescope at 5.0, 8.4 and 15.3 GHz on 10 Sept 2005, for a total of 24 hours. The 
three radio galaxies were each observed for Ri2.5 hours at each frequency. 3C84 was used as 
the instrumental polarization calibrator, assuming it to be unpolarized at these frequencies. 
The compact AGN 1156+295 was used as the electric vector position angle (EVPA or x) 
calibrator. Five to six several-minute scans of each of these sources at each frequency were 
interspersed throughout the observing schedule. 

The data were reduced using standard procedures in the Astronomical Image Processing 
System (AIPS) software package, using Los Alamos as the reference antenna at all stages 
of the calibration. Amplitude calibration was performed using the system temperatures 
and gain curve provided and the AIPS task APCAL, while the group delay and phase-rate 
calibration were performed using the AIPS task FRING. The cross polarized delays were 
determined using the procedure VLBACPOL. Finally the residual delays and delay rates 
were derived using FRING. The instrumental polarization leakage terms were determined 
using the final self-calibrated map of 3C84 in the task LPCAL. These were typically of the 
order of 1 — 2% at 5.0 and 8.4 GHz and 1 — 3% at 15.3 GHz for all the antennas. 

We determined the absolute EVPA calibration by comparing the total VLBI-scale po- 
larization of 1156+295 with integrated polarization obtained from monitoring databases 
maintained by the NRAO and the University of Michigan Radio Astronomy Observatory 



http: / /www. via. nrao.edu/astro/calib/polar/ 
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(the latter data kindly provided by H. and M. AUer). Unfortunately, there were no inte- 
grated polarization measurements with good signal-to-noise ratio available on the actual 
date of our VLBI experiment, making it necessary to interpolate these measurements to our 
VLBA date. This was straightforward for the 5 GHz EVPA calibration, since the NRAO 
database contained 5 GHz measurements on 4 and 26 September that virtually exactly co- 
incided (differing by only 1° in the polarization angle %, with errors in x being 0.5 — 1.5°). 
Accordingly, we adopted the integrated 5 GHz EVPA for 1156+295 to be -86° = +94°, and 
consider the overall uncertainty in our 5 GHz EVPA calibration to be no larger than 1 — 2°. 

NRAO measurements were also available at 8.4 GHz on these two dates, but showed a 
larger difference: x = 81° ± 1° on 4 September, and 56° ± 3° on 26 September. In addition, 
the UMRAO database had an 8 GHz measurement on 24 September with x — 75° ± 8°; this 
latter measurement is of limited use due to its large uncertainty, although it suggests that 
the EVPA of 1156+295 was not wildly variable near the time of our VLBA observations. If 
the 25° change between the integrated EVPAs for the two NRAO measurements was fairly 
smooth, the 8.4 GHz EVPA on our observing date would have been 74°. and we adopted this 
as the integrated 8.4 GHz EVPA. This value is obviously somewhat uncertain; we accordingly 
consider the overall uncertainty in our 8.4 GHz EVPA calibration to be 3 — 4°. 

The NRAO polarization monitoring database does not include measurements at or near 
15.3 GHz; the UMRAO database had two measurements near our VLBA observing date, but 
corresponding to very different x values, each with relatively large uncertainties. Therefore, 
we explored two different ways to obtain a estimate of the integrated 15.3 GHz x for our 
VLBA date: (i) obtain a linear fit of polarization angle x versus observing wavelength for 
the 4 September 2005 NRAO data at 5, 8.4, and 22 GHz and use the resulting i?M together 
with the reliably inferred 5 GHz integrated x for 1156+295 on 10 September to derive a 
15.3 GHz X value for that date; (ii) interpolate the NRAO 22 GHz x values for 1156+295 
for 4 September (63°) and 26 September (99°) in the same way that we did for the 8.4 GHz 
X\ derive a x versus plot using the interpolated 5, 8.4 and 22 GHz polarization angles 
for 10 September, and use the resulting RM for 10 September together with the reliably 
inferred 5 GHz x value for 1156+295 for 10 September to derive a 15.3 GHz x value for that 
date. In both cases, the plots of x versus A^ could be approximated by linear fits within the 
errors. Approach (i) yielded X15.3 ~ 68°, while approach (ii) yielded X15.3 ~ 72°. Realizing 
that some uncertainty in our EVPA calibration at 15 GHz was inevitable given the limited 
calibration information, but encouraged by the fact that these two estimates were similar, 
we adopted the average value of 70°. We estimate the overall uncertainty in our 15.3 GHz 
EVPA cahbration to be 5 - 6°. 

For reference, a total EVPA uncertainty of 6° between 5 and 8 GHz translates into an 
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uncertainty in the corresponding RM of about 45 radm~^. It is important to note that the 
apphcation of an incorrect EVPA cahbration will add the same offset to the observed EVPAs 
at every pixel in the map in which polarization is detected. Therefore, this can certainly bias 
the derived rotation measures, but cannot give rise to spurious i?M gradients. Thus, although 
the absolute rotation measures that are derived will have some added uncertainty due to the 
uncertainty in our EVPA calibration, this is not propagated to an added uncertainty in the 
RM gradients. 

The AIPS tasks IMAGR and CALIB were used iteratively for the imaging and self- 
calibration, respectively. Stokes Q and U maps were also made using IMAGR and the final, 
fully self-calibrated data. Pairs of the final Stokes I images at neighboring frequencies (5 
+ 8 GHz and 8 + 15 GHz) were made with identical beams, and used to derive spectral- 
index maps using the task COMB. Since information about the absolute position of the 
VLBI structure is inevitably lost during the phase self-calibration, it is necessary to first 
align the two / maps. The hybrid-mapping processes essentially aligns the VLBI cores, 
which is not th e correct absolute alignment, since the core position is frequency-dependent 



(lLobanovlll998l ). We determined the relative alignment between pairs of / images based on a 



compa rison of their optically thin regions of emission using the program of lCroke fc Gabuzda 



( 120081 ): note that the algorithm used does not rely on the presence of distinct components 



at each frequency, and so is well suited to the images considered here. 

The Stokes Q and U maps were used in the task COMB to obtain maps of the polariza- 
tion intensity and polarization position angle. These polarization-angle maps at the different 
frequencies were used in either COMB or the task 'RM' to obtain two-frequency or three- 
frequency rotation-measure maps, respectively. Prior to obtaining the i?M maps, versions of 
the Q, U, and polarization-angle maps were made with identical beams. Polarization values 
with S/N less than 3 were blanked in COMB to produce the polarized intensity maps. The 
polarization angle maps were restricted (by blanking in COMB) to have output errors less 
than ^ 10° . The three-frequency i?M map of 3C78 was made using the AIPS task 'RM', with 
unblanked polarization angle maps and a^M blanking in 'RM' such that aRM < 55 rad m~^. 
Estimates for the peak surface brightness and the RMS noise in the final total and polarized 
intensity maps were obtained using AIPS procedures TVMAXF and TVWIN & IMSTAT, 
respectively. Depolarization maps were made using the degree of polarization maps at 5, 8 
and 15 GHz for 3C78 and 5 and 8 GHz for 3C264. T he polarization intensity maps were 



first corrected for Ricean bias (e.g.. iKilleen et al.lll986l ) by using the AIPS task POLCO. 
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Results 



The total intensity radio images with polarization electric vectors superimposed are 
presented in Figs. [THSl the 5—8 GHz spectral-index maps in Fig. H] and rotation- measure 
maps and plots in Figs. [SHHl Table [1] lists some parameters pertaining to the final images. 
We discuss our results for each radio galaxy below. 

3C66B: Our observations reveal tentative polarizat ion close to the rad io core only at 



8 GHz, consistent with the epoch-2002 o bservations of iKharb et a. 



detect any polarization at 5 and 15 GHz. ISimard-Normandin et al 



(12005 ). We failed to 



(jl981bl ) have listed an 



integrated rotation measure of —67 ± 3 rad m~^, corresponding to a rotation of about 5° at 
8 GHz. The total intensity images at 5 and 8 GHz show a structure resembling a counterjet, 
as previously reported by iGiovannini et al.l (120011 ): indeed, our spectral- index map shows 
the core region itself to be optically thick, but the emission to the southwest of the core 
to be optically thin. Using the jet to counterjet surface brightness ratio, Rj = {ji^-^-^y 
(Appendix A, lUrry fc Padovanil 119951 ) at a distance of ~2 mas from the core, we obtain 
P cos 6'=0.47 for Rj ~ 22 and p{= 2 + a) = 3. We can also estimate the doppler factor, 6, 
using the ratio of the observed radio core prominence parame ter, R^., to the intrinsic radio 



core prominence, i?*"*, {6 = (-^i)^ , lUrry fc Padovanil 1 19951 ). keeping in mind the caveat 



that Rc is a statistical measure of orientation. Assuming i?*" to be 0.33 and i?c=0.05 for 



see 



3C66B 
of - 60° for 3C66B. 



Kharb fc Shastrill2004l ). the above derived f3 cos 9 gives a jet inclination angle 



3C78: Among the three FRI radio galaxies considered here, 3C78 shows the most exten- 
sive polarization along its radio jet. We detect polarization in the jet at ~3 and 6 mas from 
the VLBI core at 5, 8 an d 15 GHz. The 8 GHz image looks very similar to those from our 
epoch-2002 observations ( Kharb et al.lboosi ) with the exception of a polarized core. We did 
not detect core polarization at any frequency in our present data. This is likely to be due to 
variability in the core, since our present image sensitivity is an improvement by a factor of 
~1.5, over the 2002 observations. The total intensity images at 5 and 8 GHz show evidence 
for a counterjet, and this appears to be borne out by our spectral-index map, which shows 
the core to be optically thick, but the emission to the southwest of the core to be optically 
thin. The jet to counterjet ratio, Rj varies from ~ 32 at ~ 1 mas to 



■^22 at ~2.5 mas. This 
=0.33 and Rc=0.38 for 

3C78 (see iKharb fc Shastrill2004l ). the derived f3 cos 9 results in a jet inclination angle of 



results in (3 cos 9 varying from .52 to 0.47. Again assuming R[ 



int. 
c 



~ 50° for 3C78. 

3C78 is the only one of the three radio galaxies for which the rotation measure could be 
estimated using all the three frequencies. The rotation measure in 3C78 ranges from about 
-1-200 rad m~^ to —400 rad m~^. We observe an i?M difference of > 300 rad m~^ transverse 
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to the inner VLBI jet over a region that is ~5 mas (two beam-sizes) across (Fig. [5]). As the 
EVPA rotates through almost 90° at the lower edge of the jet (Fig. [5]), the low energy thermal 
electrons that produ ce the rota t ion cannot be comp letely mixed in with the synchrotron- 
emitting jet volume (jBurnlll966l : ICiofE &: Joned 119801 ). and are therefore external to it. We 
discuss this further in §5. Figure [9] shows the polarization map of 3C78 after correcting for 
the RM. Since these regions are optically thin, the jet magnetic field is inferred to be aligned 
to the jet direction, in the region of the i?M gradient. However, the magnetic field becomes 
orthogonal to the jet flow in the next jet component, and appears to remain orthogonal 
substantially further from the core (Fig. [2|, top left). 

3C264: Our observations reveal polarized emission at 5 and 8 G Hz in two regions i n 
the jet ~1 and 4 mas away from the core, similar to that observed by iKharb et al.l (120051 ). 
We could obtain an i?M map using only the 5 and 8 GHz maps. Therefore, the polarization 
angles are subject to ±n7r ambituities, which are impossible to resolve. The two-frequency 
RM across the inner-jet varies from ~ -|-250 rad m~^ to ~ —300 rad m~^. Note that 
an ambiguity of ±180° between 5 and 8 GHz could result in an RM change of > ±2500 
rad m~^. However, we find that the two-frequency RM in 3C264 is of the same order as 
the three-frequency i?M in 3C78 and other radio galaxy jets (Table [2]). Figure [H] shows the 
polarization map of 3C264 with the effects of the tentative two-frequency RM removed. If 
the two-frequency i?M values are correct, the jet magnetic field is roughly aligned with the 
VLBI jet direction. 



4.1. Parsec-scale RM in different AGN classes 



Prior to this study, parsec-scale rotation me asures had been obtained for only the 
broad-line FRII radio gal axies 3C111 and 3 C120 (jZavala fc Taylorl |2002| ). the narrow- line 



FRII radio galaxy 3C166 (jTavlor et al.ll200lf) and the FRI radio galaxies M87 (jJunor et al. 



200ll : IZavala fc Tavlorl 120021 ) and 3C84 flTavlor et al.l 120061 ). It is worth pointing out that 
3C120 is a peculiar radio source. Although it has sometimes been classified as an FRI (e.g., 
Owen fc Laing|ll989l) . it possesses a quasar-like prominent broad-hne spectrum ( jTadhunter et al, 
19931 ). Akin to FRII sources, it lacks co unter iet emission a nd exhibits superluminal jet mo- 
tion on scales of a few hundred parsecs (j Walker et al.lll987l ). Table [2] lists all available VLBI 
RM estimates for the radio galaxies. A non-uniform pars ec-scale rotation mea sure distribu- 



tion h as also been observed in 3C166 (ITaylor et al.ll200ll ). A recent paper by iGomez et al. 
( 120081 ) demonstrates a clear i?M gradient across the jet in 3C120. 



The absolute rotation measure estimates for the parsec-scale jets of quasars and BL Lacs 
lie typically in the range of 100—500 rad m~^ ( IZavala fc Taylorl 120041 ). This is interesting 
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Table 1. Observed Parameters 



Source 


z 


Scale 
pc/mas 


Preq 
GHz 


^peak 

mJy/bm 


hot 
mJy 


hms 

/xJy/bm 


Ppeak 

/xJy/bm 


p 

^ rms 

/xJy/bm 


3C66B 


0.0212 


0.424 


5.0 


113.9 


186.2 


43.4 




34.3 








8.4 


110.1 


193.2 


46.9 


636 


37.5 








15.3 


97.4 


160.9 


87.8 




66.6 


3C78 


0.0286 


0.567 


5.0 


342.1 


547.8 


70.4 


846 


40.4 








8.4 


330.9 


548.1 


64.7 


993 


39.8 








15.3 


242.7 


416.5 


183.7 


720 


97.1 


3C264 


0.0217 


0.433 


5.0 


144.4 


191.6 


47.5 


578 


39.8 








8.4 


139.7 


186.8 


51.8 


364 


39.4 








15.3 


106.3 


142.1 


96.3 




70.2 



Note. — Col.l: Source name. Col. 2: Redshift. Col.3: Parsecs corresponding to 1 
mas in source. Col.4: Observing frequency in GHz. Col.5: Peak surface brightness 
in total intensity (/) map. Col. 6: Total VLBI flux density. Col.7: RMS noise in / 
map. Col. 8: Peak surface brightness in polarized intensity (P) map. Col. 9: RMS 
noise in P map. 



- 10 - 



Table 2. Radio galaxies with measured parsec-scale rotation measures 



Source 


FR 


z 


Int. RM 


RM 


Comment 


Ref 


Name 


Type 




rad m~^ 


rad m~^ 






3C78 




0.0286 


+14 


-400, +200 


in jet at r= 1.7, 4.3 pc 


1 


3C264" 




0.0217 


+16 


-250, +150 


in jet at r= 2.3, 1.9 pc 


1 


3C84 




0.0175 


b 


+7000 


in jet at r= 5 pc 


2 


M87 




0.0043 


+872 


-4500 


in jet at r= 2 pc 


3 


3C111 


11 


0.0485 


-19 


-750, -200 


in jet at r= 2.8, 4.7 pc 


3 


3C120 


11 


0.0330 


-3 


+ 100^ 


in jet at r= 1 pc 


3 










+1500*, +5000 


in jet at r= 1, 2 pc 


4 


3C166 


11 


0.2449 


+ior 


-2300*, -310* 


in jet at r < 3.8, 7.6 pc 


5 



Note. — Comparison of FRl and FRll radio galaxies with parsec-scale RM mea- 
surements available at present. " Parsec-scale RM was obtained using only two fre- 



Simard-Normandin et al.l (Il981bl ). * Unpolarized in 



Simard-Normandin et al. 


fl981a) 


Simard-Normandin et al. 


f 1981a) 



Integrated RM was estimated from values listed in 
Col. 5 lists the average i?M at a pro jected distance r fro m the core. | Core RM of +2080 
rad m^^ was detected in 3C120. | iGomez et al.l (120081 ) have suggested a time variabil- 
ity in the jet RM of 3C120. * Rest-frame rotation measure, obtained b y multiplying 



RM b y (1 + 2:)^=1.55 for 3C166 . Ref erences: 1 - Presen t pap er. 2 - iTavlor et al 
J2006h 3 - IZavala fc Tavloil J2002h . 4 - IComez et aP J2OO8I 1 5 - iTavlor et all J200lh . 
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because the absolute parsec-scale jet RM values of the FRI and FRII radio galaxies studied 
so far, although few, fall in the same range. Further, the galaxies studied so far seem to reside 
in widely varying environments, from the centre of clusters (e.g., 3C264) to being relatively 
isolated (e.g., 3C78). Therefore the similarity in the i?M values across the parsec-scale jets 
of different classes of AGN, suggests that the Faraday rotating medium is perhaps intrinsic 
and closely associated with the jet itself. 



5. The Faraday rotating medium 
5.1. Galaxy ISM 



I t has been clear frora the earliest VLBI Faraday rotation measu rements (e.g.. iTaylor 



19981 : iNan et al.lll999l . l2000l : iGabuzda et al.ll200ll : iRevnolds et al.ll200lh that the parsec-scale 
i?M distributions are usually non-uniform, indicating that there must be a substantial compo- 
nent of Faraday rotation that is intrinsic to the AGNs , rather than arising in our own Galax y 
or on larger scales in the host galaxy of the AGN (e.g.. lTaylorlll998l : IUdomprasert et al.lll997l ). 
The RM varies on scales of a few parsec in 3C78 and 3C264. The RM gradient across the 
VLBI jet of 3C78 is observed over a ~ 5 mas (~ 3 pc) region. Therefore it is difficult to assign 
the role of the Faraday rotating medium to the IS M in our Galaxy or the radio host galaxy. 
Recent X-ray observations of elhptical galaxies (see lMathews &: Brighentill2003l ) suggest that 
the electron density of the hot gas in the centers is typically ~ 0.1 cm~^ and falls off with 
radius (r) as oc 7--i-25±o.25_ Assuming that the RM is produced in the hot gas prevalent 
on scales of a 100 pc, it would produce an RM >40,000 rad m~^ for a uniform magnetic 
field, where the line of sight magnetic field strength is equal to the 'equipartition' magnetic 
field (say ~5 mG, see §5.2). This is clearly not observed. However, if the magnetic field was 
not uniform and there were 'N' reversals in the field direction along the line of sight, then 
the observed RM of ~200 rad m~^ would require N> 40,000. Such a large number of field 
reversals would however result in a net depolarization on parsec-scales, which does not seem 
to be the discussed below. 



5.2. Radio Jet and Surrounding Medium 



Considering the medium surrounding the parsec-scale jet : the broad- line region (BLR) 
has been inferred to have extents typically less than a parsec ( lUrry fc Padovanilll995l ). BLR 
clouds could therefore be ruled out as a potential Faraday rotating medium for the jet region 
lying a couple of parsecs away from the core. The narrow-line region (NLR) could however, 
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extend to hundreds of parsec s and have a large covering factor (~ 30 — 50%, iNetzer fc Laor 
19931 : iRowan-RobinsonI Il995l ). The NLR clouds could therefore in principle contribute to 
the observed rotation measures. However, the short timescale variability observed in the 
parsec-scale RM suggests that NLR clouds or the intercloud gas cannot be the primary 
Farad ay rotating mediurn, but rather this med ium is closely linked with the relativistic jet 
itself ksada et aliboosl : IZavala fc Tavto^boosl ). 



We created depolarization maps for 3C78 and 3C264 using the relation DP = where 
mi, rrih are the degree of polarization maps at the lower and higher frequencies, respectively. 
We made DP maps for 3C78 between 5 and 8 GHz (polarized intensity maps were blanked 
below lOcr) and between 8 and 15 GHz (polarized intensity maps were blanked below 3a). 
The mean 5-8 GHz and 8-15 GHz depolarization values in 3C78 in the region ~ 3 mas away 
from core were DPg = 0.98 ± 0.11 and DPfr, = 0.79 ± 0.15, while in the region ~ 6 mas 
away they were DP^ = 0.60 ± 0.07 and DPf^^ = 0.64 ± 0.09. We note that the DPf^, value 
in the region closer to the core (~ 3 mas) was obtained over a much smaller region than 
DP^. Overall these large DP values suggest that there was little depolarization at longer 
wavelengths in 3C78. The depolarization map for 3C264 was made using the degree of 
polarization maps at 5 and 8 GHz (polarized intensity maps were blanked below 5a). The 
mean depolarization was DP^ = 0.44 ± 0.11 around ~ 1 mas away from the core. 

We noted in §4 that in 3C78, the polarization angle seems to rotate by almost 90° at the 
lower edge of the jet (Fig. which indicates that the Faraday rotating medium is not mixed 
with the main synchrotron-emitting jet. Furthermore, the low depolarization suggested by 
the DP estimates and the strong dependence of polarization angles on the square of the 
observing wavelength, supports the picture of a Faraday screen with a relatively constant 
Faraday depth. Since the hot gas observed in X-rays in fact requires a highly tangled 
magnetic field to produce the observed RM, which does not seem to be borne out by the 
depolarization measurements, we suggest that an optically thin layer surrounding the jet 
that has sufficient thermal electrons mixed in with any synchrotron-emitting electrons could 
be a good candidate for the Faraday screen. 

Sheath s around jets have been suggested to be present in both FRI and F RH jets on both 



kiloparsec- (IHardcastle et al.lll996l : iKatz-Stone et al.lll999l ) and parsec-scales (lAttridge et al. 



1999 : Kharb et al. 20081 ). Either shearing due to the jet interaction with the ambient medium 
or the presence of an intrinsic to roidal or hel i cal magnetic field cou ld result in these sheaths 
having ordered magnetic fields (jLaingj Il98ll : iLyutikov et al.l 120051 ) . Although a sheath re- 
sultin g from jet entrainment has often been proposed for FRI jets ( iDe Younglll984l : iBicknell 
19961 ). sheaths co uld also result from t he diffusion of particles out of the beam into the 
external medium (jO'Dea &: OwenI 119871 ). or from material associated with the inner part 
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of a disk wind (iBlandford fc Levinsoru Il995l : iHanasz fc Soil Il996l ) . The fact that we observe 
similar parsec-scale jet rotat ion measures in FRIs and FRIIs, where the levels of entrainment 
are expected to differ (e.g., iDe Younglll993l ). appears to provide support to the latter two 



models, under the assumption that jet-like magnetic field lines also thread the sheath region. 

The most natural interpretation fo r the RM gradient observed in 3C78 is that it is due to 
a helical or toroidal jet magnetic field (IBlandf or dl 1 1 99 3l : lLaing]|l996l ). This is also consistent 
with the increase in the degree of polarization towards the edge of the jet (FigE]). The 
asymmetry in the i?M values is likely to be the combined e ffect of the jet orientation and the 
pitch angle of the field (for example, see lAsada et al.ll2002l ). Similar RM g radients have been 



observed in a number of other AG Ns, including several BL Lac o bjects ( lAsada et al.l 12002 



Gabuzda et al.ll2004l : IZavala fc Tavloii l2005l : iGabuzda et al.ll2008[ ). Thus, the i?M gradient 
together with the large rotation in the polarization angles and the low depolarization supports 
the idea that the Faraday rotating medium is a 'sheath' around the parsec-scale jet, which 
is threaded by a helical magnetic field that is intrinsically associated with the jet. 

The electron density in the medium producing the Faraday rotation can be derived 
under the assumption t hat the condit ion of 'equipartition' betwee n magnetic field and 
electron energv density (jBurbidgd Il959l ) holds in the radio jet (e.g., iFeigelson et al.l Il995 



Harris fc Krawczynskill200^ . Using the 5 GHz jet flux density and volume we obtained the 



magnetic field streng th, electron energy de nsity and pressure (Table [3]) under the 'equipar- 
tition' condition (see lO'Dea fc Owerull987l ). For these calculations, we assumed the ratio of 
the relativistic proton to relativistic electron energy to be unity. The radio spectrum was 
assumed to extend from 10 MHz to 100 GHz. An average spectral i ndex of 0.8 and a volume 
filling factor of unity was adopted for the jet. The constant C12 (jPacholczykI Il970l ) which 
depends on the spectral index and frequency cutoffs was taken to be 6.5 x 10^. 

Assuming that the Faraday rotating medium was external to the synchrotron-emitting 
jet volume but close to it, with the magnetic field strength being similar to the jet 'equipar- 
tition' 5-field of ~ 5 mG a canonical rest-frame rotation measure of 200 rad m~^ and a path 
length of 1 pc, we derive a lower limit to the electron density which is a few times 0.01 cm~^ 
for the medium in 3C78 and 3C264. If we assume that the path length is in fact the width 
of an outer sheath layer, which is (say) about 10% of the jet width (^5 pc, estimated from 
the 5 GHz maps), we obtain an electron density of ~ 0.1 cm~^, which is of the same order 
as the electron density of the hot gas in the galaxy centers. Therefore if the sheath was 
a 'mixing' layer between the jet and the surrounding hot gas, it could certainly have the 
required electron density, magnetic field strength and path length, to produce the observed 
RM. 



Although helical jet S-fields provide an attractive explanation for the transverse RM 
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gradient and increase in fractional polarization at the jet edges in 3C78, alternative explana- 
tions, such as an inhomogeneous surrounding medium (e.g., accretion disk wind) that gives 
rise to a local enhancement in the RM on one side of the jet, cannot be completely ruled 
out. The absence of a clear transverse RM gradient in 3C78 further from the core may 
provide evidence for an inhomogeneous ambient medium, whether or not this jet carries a 
helical S-field. The jet of 3C264 does not show any clear gradients in our two-frequency 
RM map, suggesting that either this jet does not carry a helical S-field, or that, if the jet 
does carry a helical field, this field is distorted by interaction with the surrounding medium. 
Higher-resolution low-frequency observations that are sensitive to Faraday rotation but bet- 
ter resolve the jet in the transverse direction, such as may be provided by future space VLBI 
projects, may help resolve these questions. 



5.2.1. Clues from optical observations 



The three FRIs under study have been observed with optical imaging polarimetry with 
the HST at a resolution of ~ 0.1" (tr anslating to scales of h undreds of par s ecs fo r these 



sources) by iPerlman et al.l (120061 ) and ICapetti et al.l (120071 ). iPerlman et al.l (120061 ) found 
that 3C66B, 3C78 and 3C264 showed the simplest magnetic field structures among the six 
nearby radio galaxies studied. Interestingly, both 3C78 and 3C264 (sources with significant 
parsec-scale polari zation) exhibit an in crease in the degree of polarization along the edges 
of the optical jet. IPerlman et al.l (120061 ) suggest that such a magnetic field structure could 
arise due to a spine-sheath jet structure, with the sheath forming due to the interaction of 
the jet with the surrounding medium, which stretch the magnetic field lines due to shear. 
We note that a spine-sheath st ructure could also sim ply result from a helical jet magnetic 
field, with or without shearing (ILyutikov et al.ll2005l ). 



A 'potential' problem with the jet interaction picture could be that on going down 
from ~ 0.1" to mas-scales, we detect polarization along the centre of the jet (in 3C78 and 
3C264) but not at the edges as expected from a 'sheath'. This could either be due to 
insignificant interaction/entrainment on parsec-scales, or, the existence of a helical magnetic 
field extending from parsecs (where its presence is revealed by the observed transverse RM 
gradient) to hundreds of parsecs (where its presence is revealed by the observed transverse 
optical polarization structure). This helical field could originate due to the rotation of the 
accretion disk plus jet outflow, and propagat e outward with the jet throughout its length. 
Another possibility is the scenario proposed by lAsada et al.l (120081 ) for the quasar NRAO140, 
viz.., that the jet has a spine-sheath structure with helical magnetic fields threading both the 
jet components, but the field is tightly wound in the sheath (resulting in the RM gradient) 
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and loosely wound in the spine (leading to aligned magnetic fields in the inner jet). This 
scenario is also consistent with the findings in 3C120 and 3C166. Interestingly, such a 
magnetic field geometry has been known to result in the matter + Poynting flux domin ated 
jets produced in the 3-D general relativistic MHD simulations of iDe Villiers et al.l (120051 ) and 
Hawley fc KrolikI (120061 ) . Such a magnetic field geometry could result in net depolarization 
at the jet edges which have transverse fields as opposed to the dominant spine which have 
aligned fields, thereby resolving the 'potential' problem. Significant interaction/entrainment 
on scales of a few hundred parsecs, along with magnetic fiux conservation could reproduce 
the tranvers e -B-field spine and l ongitudinal i?-fiel d sheath observed in kiloparsec-scale radio 
galaxy jets (ICanvin et al.ll2005l : iLaing et al.ll2006l ). 



Another interesting feature that emerges from optical observations is the presence of 
optical jets in the majority of radio galaxies with detected parsec-scale rotation measures, 
i.e., detected parsec-scale polarization over two or more radio frequencies, which in turn 
could be the result of low depolarization between frequencies. With the exception of 3C111 
and 3C166, all the radio galaxies listed in Table 2 h ave optical jets. Ho wever, even 3C166 
exhibits an optical 'tail' along the radio jet direction (Ide Koff et al.lll996l ). Both 3C111 and 
3C166 also exhibit u nresolved iJ^T nuclei which have been suggested to be the unresolved 
bases of optical jets (IChiaberge et al.ll2000l ). A possible connection between the presence of 
optical emission in the jet and detectable parsec-scale polarization in the radio needs to be 
examined further with a greater number of sources. 

However if there was such a connection, it can be easily understood. Optical jet emission 



Sparks et al. 
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(.e.g. 

emission, viz.., 3C111 and 3C120, are broad-line radio galaxies, and are likely to exhibit 
Doppler boosting in their jets due to being oriented at relatively small angles to line of sight. 
Furthermore, the kpc-scale radio core prominence values of the three FRIs under study 
indicate that the most polarized source, 3C78, is also the most core-dominant, while the 
least polarized source, 3C66B, is the least core-dominant. This underscores the importance 
of Doppler boosting in these parsec-scale jets. Particle re-acceleration could occur in the jet 
shear layer or sheath, as the jet interacts with the surrounding medium, thereby producing 
a net deceleration and converting bulk kinet ic energy into high energy optical radiation 
(lOwen et al.lll989l : IStawarz fc Ostrowskil 120031 1. 
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6. Summary and Conclusions 

We have observed three FRI radio galaxies, viz., 3C66B, 3C78 and 3C264, with polar- 
ization sensitive VLBI at 5, 8 and 15 GHz. We detected polarization in all three and obtained 
rotation measures across the parsec-scale jets of 3C78 and 3C264. RM was obtained using all 
the three frequencies in 3C78 but only two frequencies in 3C264, which is therefore subject 
to uncertainty due to ±n7r ambiguities that cannot be resolved. To summarise, 



1. Polarization was detected in the jets of 3C66B at 8 GHz, in 3G78 at all the three 
frequencies and in 3G264 at 5 and 8 GHz. 

2. The total intensity maps confirm the detection of a counterjet in 3C66B and indicate 
the presence of one in 3C78. The 5-8 GHz spectral index maps reveal an optically 
thick core and an optically thin jet in all sources. The counterjet emission in 3G66B 
and 3G78 is also optically thin. 

3. The parsec-scale jet rotation measure in 3C78 ranges from ~+200 rad m^^ to ~ —400 
rad m~^. The core appears to be completely depolarized. A rotation-measure gradient 
is observed across the jet in 3G78 in a region that is over two beam-sizes across. Fur- 
thermore, the degree of polarization increases along the edge of the jet. This strongly 
supports the idea of a helical jet magnetic field. 

4. The polarization angles rotate by almost 90° between the three frequencies in 3C78. 
The depolarization parameter is close to unity in the jet, suggesting low depolarization 
at longer wavelengths. This suggests that the Faraday rotation is occuring in a layer 
containing thermal plasma that is external to the main body of the jet. 

5. After correcting for the Faraday rotation, the magnetic field in 3C78 seems to be 
aligned with the jet direction close to the core, but becomes orthogonal further down 
the jet. 

6. The two-frequency parsec-scale rotation measure in the 3C264 jet varies from ~-|-250 
rad m~^ to ~ —300 rad m~^. The core is completely depolarized. 

7. The parsec-scale rotation measures in 3G78 and 3G264 are similar to those observed in 
other nearby FRI and FRII radio galaxies. Parsec-scale RM of a few hundred rad m~^ 
have also been inferred in the jets of the Doppler-beamed BL Lacs and quasars. The 
similar parsec-scale rotation measures in these AGN classes suggests an intrinsic origin, 
perhaps in the jet medium itself. 
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Based on the transverse i?M gradient, increase in the degree of polarization towards the 
edge of the jet, large rotation of the polarization angles due to the Faraday rotation, and 
the small depolarization in 3C78, we argue that a layer surrounding the jet and carrying 
a helical magnetic field forms the Faraday scree n. This sugg e stion is supported by HST 



optical polarization images. The recent work by iGomez et al.l ( 120081 ) on 3C120 is likewise 
fully consistent with our conclusions. 

The presence of optical jets in the majority of radio galaxies with detected parsec-scale 
polarization and the spine-sheath polarization structure observed in the HST polarization 
images of some of these sources suggests a simple connection: a sheath layer around the jet, 
produced either by jet-medium interaction or a helical magnetic field, is producing both. 
Jet-ISM entrainment and interaction could cause (1) mixing of thermal gas with the outer 
sheath layer of the jet which produces the rotation measure, and (2) particle acceleration 
which produces the optical jet. This connection needs to be further explored based on a 
larger number of sources. 
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Fig. 1.— Total intensity maps of 3C66B at 5.0 GHz (top left), 8.4 GHz (top right) and 15.3 
GHz (bottom) with polarization vectors superimposed. The contours are in percentage of the 
peak surface brightness and increase in steps of 2. The lowest contour and peak brightness 
for 5, 8 and 15 GHz are ±0.085% of 114 mJy beam'^ ±0.085% of 111 mJy beam'^ and 
±0.085% of 97 mJy beam~^, respectively, x vectors for 8 GHz: 1 mas = 0.3 mJy beam"-*^. 
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Fig. 2. — Total intensity maps of 3C78 at 5.0 GHz (top left), 8.4 GHz (top right) and 15.3 
GHz (bottom) with polarization vectors superimposed. The contours are in percentage of the 
peak surface brightness and increase in steps of 2. The lowest contour and peak brightness 
for 5, 8 and 15 GHz are ±0.042% of 342 mJy beam-\ ±0.042% of 331 mJy beam-^ and 
±0.17% of 242 mJybeam"^, respectively, x vectors for 5, 8 and 15 GHz: 1 mas = 0.25, 0.3, 
and 0.8 mJy beam~^, respectively. 
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Fig. 3.— Total intensity maps of 3C264 at 5.0 GHz (top left), 8.4 GHz (top right) and 15.3 
GHz (bottom) with polarization vectors superimposed. The contours are in percentage of the 

peak surface brightness and increase in steps of 2. The lowest contour and peak brightness 
for 5, 8 and 15 GHz are ±0.085% of 144 mJy beam"\ ±0.085% of 140 mJy beam"\ and 
±0.17% of 106 mJy beam^^, respectively, x vectors for 5 and 8 GHz: 1 mas = 0.25 and 0.3 
mJy beam~^, respectively. 
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Fig. 4. — The 5—8 GHz spectral index map of (top left) 3C66B made with a 1.5 mas circular 
beam, (top right) 3C78 and (bottom) 3C264, each made with a 2 mas circular beam. The 
5 GHz total intensity contours are superimposed. Note the optically thin counterjet emission 
in 3C66B and 3C78. 
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Fig. 5. — The 3-frequency rotation measure image of 3C78 with 5 GHz total intensity 
contours superimposed. An area covering four pixels (pixel size=0.2 mas) was used to create 
the inset x-^^ plots. The RM values in the inset plots and colour scale are in units of 
rad m~^, with aRM < 55 rad m~^. Rotation measures in the rest-frame of the source are 
larger by a factor (1 + z)'^ or 1.058. 
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Fig. 6. — 5 GHz total intensity contours with degree of polarization vectors superimposed for 
3C78 (left) and 3C264 (right). 1 mas ticks correspond to 8% and 50% for 3C78 and 3C264, 
respectively. Degree of polarization clearly increases towards the southern and eastern edge 
of the 3C78 and 3C264 jets, respectively. 

Table 3. Derived Parameters 



Source 


Lrad 


E ■ 

-'-'mm 




P . 

^ mm 


Tie 




erg/s 


erg 


mG 


dynes/cm^ 


cm~^ 


3C66B 


3.9E+40 


6.9E+51 


8.1 


6.1E-06 




3C78 


2.1E+41 


4.6E+52 


7.1 


4.7E-06 


0.03 


3C261 


1.3E+10 


1.2E+52 


5.7 


3.0E-06 


0.01 



Note. — Col.l: Source name. Col. 2: The averaged total 
radio luminosity. Col. 3 & 4: particle (electrons and protons) 
energy and magnetic field strength at the minimum pressure 
Pmin, obtained assuming the 'equipartition' condition. Col. 5: 
Minimum Pressure. Col. 6: Electron density derived for an 
RM of 200 rad m~^ using the 'equipartition' 5-field strength 
and path length dl—1 pc. 
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Fig. 7. — The RM gradient across the jet in 3C78. The left panel zooms in and displays 
a narrower range of RM where the RM gradient is clearly visible perpendicular to the jet 
close to the core. The RM changes from ~ —450 to ~ —50 rad m~^ over a region 5 mas 
(two beam-sizes) across the jet (Inset A). Inset B displays the rotation measure along the jet 
direction. Rotation measures in the rest-frame of the source are larger by a factor (1 + z)'^ 
or 1.058. 
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Fig. 8. — The rotation measure image of 3C264 with 5 GHz total intensity contours su- 
perimposed. Only the 5 and 8 GHz images were used to derive the RM. The colour scale 
is in units of rad m~^ with (T/jm being typically < 50 rad m~^. Rotation measures in the 
rest-frame of the source are larger by a factor {1 + z)"^ or 1.044. 
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Fig. 9. — The polarization maps with 5 GHz total intensity contours of 3C78 (left) and 
3C264 (right) with the effects of i?M removed. Circular beams of 2.5 and 2.0 mas were used 
for the 3C78 and 3C264 maps, respectively, x vectors for 3C78 and 3C264: 1 mas = 0.6 
and 0.4 mJy beam~^, respectively. The intrinsic EVPA for the jet in 3C78 changes from 
nearly longitudinal to nearly transverse to the VLBI jet direction, while it remains roughly 
transverse throughout in 3C264. 
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